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Point  groups:  For  the  point  groups  (crystal  class)  the  standard 

Schonflies  labels  (0^,  C2h*  Dj,  etc.)  are  used. 

Labeling  of  irreducible  representations:  For  the  cubic  group  (0^) 

the  BSW  notation1  is  used.  for  the  tetragonal  group  (D^)  is 
identical  to  of  BSW.  For  the  hexagonal  group  (Dgh)  and  for  the 
additional  groups  listed  in  Table  C-l  the  notation  is  essentially  that 
of  Koster  et  al.2;  we  have  replaced  with  P^.  The  correspondence 

between  frequently  used  labeling  systems  for  the  cubic,  hexagonal,  and 
tetragonal  groups  is  shown  in  Appendix  A. 

INTRODUCTION 

A useful  check  on  the  crystal  potential  used  in  an  electronic 
bandstructure  calculation  may  be  obtained  from  a first-principles 
determination  of  deformation  potentials  if  these  can  be  checked  against 
experimentally  determined  values.  It  is  well  knownJ»  that  strains 
associated  with  a change  in  the  crystal  symmetry  provide  a particularly 
critical  test  of  the  original  potential. 

1L.P.  Bouckaert,  R.  Smoluchowski,  and  E.  Wigner,  Phys.  Rev.  50,  58  (1936) 
2"Properties  of  the  Thirty-Two  Point  Groups",  G.F.  Koster,  J.O.  Dimmock, 
R.G.  Wheeler,  and  H.  Statz  (M.I.T.  Press,  1963). 

3d.  Shoenberg  and  B.R.  Watts,  Philos.  Mag.  15^,  1275  (1967). 

4R.  Griessen  and  R.S.  Sorbello,  J.  Low  Temp.  Phys.  237  (1974). 
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The  motivation  behind  the  use  of  perturbation  theory  is  very  simple; 
one  wants  to  avoid  calculating  the  small  energy  shifts  AE  from  a 

AE  - E1  - E°  (1) 


type  of  calculation. 

In  the  unperturbed  crystal  the  energy  levels  are  determined  from  a 
solution  of  the  SchrSidinger  equation, 

C-  V2  + V (r)  ] = E°  Y^Cr).  (2) 

The  electron  wave  function  is  approximated  by 


Yn°(x)  = Ii  cj  0?(r)  (3) 

where  the  0°  are  known  functions  and  the  c?  are  determined  by  the 
variational  procedure.  The  variational  procedure  leads  to  the  set  of 
equations 


Zj  c;<0?iH°-Esi0j>"°  <4> 


where  H°  is  the  Hamiltonian  of  the  unperturbed  system  Hq  = (-V2  + V°(r)) 
and  E^  is  the  nth  energy  eigenvalue.  (We  arbitrarily  order  the  E°'s 
from  lowest  to  highest.)  In  matrix  form  the  last  equation  becomes 


= E°  S°  c , 
n n* 


(5) 


with 


so. . . <0o  | 00) 

(6a) 

H°ij  =<0°IH°|0°> 

(6b) 

As  the  0°  are  not  necessarily  orthogonal,  S®  is  not  the  identity 
matrix. 

From  the  computational  point  of  view  there  are  two  reasons  for 
utilizing  Group  Theory  in  the  E°  problem:  (1)  to  provide  automatic 

accounting  of  degeneracy  and  of  the  connectivity  between  different- 

<*  «■* 
magnitude  k points  along  the  same  symmetry  direction  in  k space.  (2) 

To  decrease  the  number  of  trial  expansion  functions  necessary  to 

obtain  good  convergence  (and  thus  decrease  computer  time) . Both  of  these 

stem  from  a general  principle  of  quantum  mechanics;  the  wave  functions 

of  a quantum  system  must  form  bases  for  irreducible  representations  of 

the  group  of  operators  which  commute  with  the  Hamiltonian  of  the  system.^ 

In  the  modified-plane-wave  formulation^'^  within  which  we  work  the 

(£)  ' s are  plane  waves  and  atomic-like  functions.  The  plane  wave  P's 

change  under  strain  due  to  containing  the  lattice  constant (s)  explicitly. 

Thus,  the  usual  first-order  perturbation  expression 

<y;iHpertiYs>  (?) 

^"Energy  Band  Theory",  J.  Callaway  (Academic  Press,  1964),  Section  1.6. 

6E.  Brown  and  J.A.  Krumhansl,  Phys.  Rev.  109,  30  (1958). 

7 

D.  Gray  and  E.  Brown,  Phys.  Rev.  160,  567  (1967). 
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i t ’I  'A  fl  . nViriiturH 


is  replaced  by 


AEn  - *ij  ci  cj  <AHij  - >5  W 

with 

Atti i 3 H.l  - H.?  (9) 

X1  13  lj  1 J 

where  J - V2  + V°(r)  + AV(r)  I 0j>-  (10) 

(See  ref  8;  in  particular  Section  II  and  Appendix  A.)  is  defined 

analagously  to  AH_.  Three  important  points  should  be  noted  here:  (1) 
AH.j  (and  AS^)  are  obtained  algebraically  by  expanding  (S^j) 
in  a Taylor  series.  (Otherwise,  one  has  merely  shifted  the  subtraction 
approach  of  eq  (1)  to  the  individual  S^  , elements.)  (2)  Both 
Hij,  S^j  and  H^9,  S^°  must  be  expressed  in  the  lower  symmetry  for  those 
cases  in  which  the  perturbation  changes  the  symmetry.  (Once  the  crystal 
is  deformed,  the  higher  symmetry  of  the  original  crystal  has  no  real 
meaning,  thus  S^j  must  be  expressed  in  the  lower  symmetry;  if  one 

then  tries  to  express  H^9,  S^9  in  the  higher  symmetry,  there  is  no  simple 
correspondence  between  the  ij  labels  of  H°,  S°  and  those  of  T?*, . ) 

(3)  To  first  order  AV(r)  equals  zero  if  both  of  the  following  conditions 
hold:  (A)  The  strain  contains  no  component  (that  strain  component  which 
preserves  the  original  symmetry);  (B)  V'r)  and  AV(r)  are  taken  as  spheri- 
cally symmetric.  (See  ref  9,  Appendix  A,  section  A.) 

^A.M.  Gray,  D.M.  Gray,  and  E.  Brown,  Phys.  Rev.  B,  11 , 1475  (1975). 

9D.M.  Gray  and  A.M.  Gray,  Phys.  Rev.  B,  14_,  669  (1976). 


For  cubic  crystals,  explicit  expressions  for  AS.^  and  , for 
hydrostatic,  tetragonal,  and  trigonal  strains  are  given  in  ref  8, 

Appendix  C. 

In  applying  this  procedure  to  specific  crystal  and  strain  types  one 
needs  to  know  what  group  is  associated  with  the  lower  symmetry.  The 
main  part  of  this  report  is  devoted  to  the  determination  of  this  lower 
symmetry  group  for  various  strains  and  to  what  simplifying  relations 
exist  due  to  symmetry  alone.  Once  the  group  associated  with  the  lower 
symmetry  is  identified,  extensive  use  may  be  made  of  existing  compat- 
ibility tables.  As  an  example,  consider  a P / level  of  a cubic 
crystal.  In  the  unperturbed  (E°)  case  ^5'  levels  are  3-fold  degenerate 
(p25/is  3-dimensional).  Now  apply  a tetragonal  strain.  Compatibility 
tables  (connecting  the  tetragonal  subgroup  to  the  full  cubic  group) 
associate  with  and  P5  of  the  tetragonal  group.  Since  P^  is  1- 

dimensional  and  P^  is  2-dimensional,  one  immediately  knows  that,  under 
this  strain,  a P^/  level  will  split  into  two  levels,  one  non-degener- 
ate  (P3)  and  one  2-fold  degenerate  (P  ).  The  unperturbed  P^  level 
will  show  up  on  P^  and  on  P^  whereas  the  different  shifts  will  show 
up  only  on  P^  orT^,  as  appropriate;  in  this  particular  example  one 
need  only  compute  one  of  these  two  shifts,  the  other  may  then  be 
deduced  immediately  from  symmetry  considerations  alone. 

SYMMETRY  THEOREMS 

The  example  just  given  illustrates  a major  simplification  due  to 
symmetry'  this  simplification  stems  from  the  following  theorem  (proved 

8A.M.  Gray,  D.M.  Gray,  and  E.  Brown,  Phys.  Rev.  B,  ]J,  1475  (1975). 

I • I 

d 
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in  ref  9) : The  Fermi  energy  Ep  and  the  sum  over  the  "star-of-k" 

of  any  given  level  belonging  to  any  given  irreducible  representation 

changes  (to  first-order)  under  strain  only  if  the  strain  contains  a 

component  (i.e.,  contains  a symmetry-preserving  component).  This 

theorem  allows  one  to  make,  from  symmetry  considerations  alone,  a number 

of  predictions  regarding  ratios  of  energy  level  shifts®.  It  should  be 

noted  that  no  more  information  can  be  obtained  from  the  Wigner-Eckart 
8 10 

theorem  * than  from  standard  compatibility  tables  combined  with  the 
"star-of-k"  theorem  above;  the  Wigner-Eckart  theorem  leads  to  the  same 
ratios  and  corroborates  the  information  that  some  shifts  are  probably 
not  equal  to  zero. 

There  are  two  important  corollaries  to  the  "star-of-k’1  theorem: 

(1)  For  cubic  crystals,  AVolume=0  is  the  necessary  (except  for 
"accidental"  cases)  and  sufficient  condition  for  AE  (etc.)=0. 

r 

(Appendix  B) . 

(2)  For  non-cubic  crystals,  AVolume=0  is  only  a necessary  condition. 
(Appendix  B). 

If  a model  is  chosen  which  constrains  the  potential  to  spherical 
symmetry  within  a muffin-tin  sphere  (after  deformation  as  well  as 
before)  the  discussion  in  ref  9 shows  that  our  theorem  applies  to 
V(r)  also;  thus,  for  such  a model,  and  for  strains  not  containing  P^, 

®A.M.  Gray,  D.M.  Gray,  and  E.  Brown,  Phys.  Rev.  B,  H_,  1475  (1975). 

^D.M.  Gray  and  A.M.  Gray,  Phys.  Rev.  B,  14_,  669  (1976). 

10"Operator  Techniques  in  Atomic  Spectroscopy,"  B.R.  Judd  (McGraw-Hill, 
1963),  pp.  39-41. 
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the  bandstructure  shifts  and  splittings  are  due  solely  to  the  change 
in  symmetry  and  to  the  explicit  change  in  lattice  constant(s),  i.e., 
there  is  a "geometric"  effect  but  no  "potential"  effect  (change  in 
bandstructure  due  to  the  change  in  V(r)  caused  by  the  change  in  lattice 
constant). 


DECOMPOSITION  OF  THE  STRAIN  AND  DETERMINATION  OF  THE  ASSOCIATED 
LOWER  SYWETRY  GROUPS 
A.  General. 

An  arbitrary  (symmetric)  strain  can  be  represented  as 


ell  e12  e13 


e12  e22  e23 


e13  e23  e33 


The  components  are  defined  as  follows:  imagine  that  three  orthogonal 

^ A ^ 

axes  fQ,  gQ,  hQ  of  unit  length  are  imbedded  securely  in  the  unstrained 

^ 

solid.  After  a small  deformation  these  axes  become  f,  g,  h with 
f = (!  + en)  fQ  ♦ e12  gQ  ♦ e13  hQ  (12) 
etc.  (This  definition  is  consistent  with  that  used  in  ref  8 and  9 
and  corresponds  to  Kittel's  notation11  except  for  a factor  i/2  in  the 


8A.M.  Gray,  D.M.  Gray,  and  E.  Brown,  Phys.  Rev.  B,  1_1,  1475  (1975). 
®D.M.  Gray  and  A.M.  Gray,  Phys.  Rev.  B,  J_4,  669  (1976). 

11" Introduction  to  Solid  State  Physics",  2nd  ed.,  C.  Kittel  (Wiley, 
1956),  Chap.  4. 
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off-diagonal  elements.)  From  eq  (11)  and  (12)  it  may  readily  be  seen 
that 


AVolume  = (Volume)  Tr  e 


(13) 


^ ^ 

by  writing  Volume  = f Xg  • h. 
We  may  write  e as 


= 7t 


(OC 


i,  oc  x oc 


(14) 


i.e.,  as  a linear  sum  of  component  strains,  each  transforming  like  one 
partner  of  one  irreducible  representation.  (In  eq  (14)  the  i index  is 
associated  with  a representation  label;  the  oc index  allows  for  cases 
where  there  is  more  than  one  independent  e associated  with  a given 
representation,  either  through  different  partners  of  a multi-dimensional 
representation  or  through  independent  functions  associated  with  a one- 
dimensional representation,  e.g.,  x^+y^  and  z^  both  associated  with 
P for  a tetragonal  crystal.)  Since  there  are  six  independent  entities 
in  eq  (11)  six  e10C  will  be  required  for  the  most  general  (symmetric) 
strain  e. 

For  any  crystal  class  (point  group)  the  irreducible  representations 
which  may  be  included  in  the  decomposition  of  an  arbitrary  strain  must 
be  contained  in  rPXyZ3s»  t*ie  symmetrized  square  of  PXyZ>  where 
PXyz  is  the  representation  (possibly  reducible)  which  has  x,  y,  and 
z as  its  basis  functions.  (A  proof  of  this  theorem  is  outlined  in  Appen- 
dix C;  Table  C-l  lists  the  irreducible  representations  contained  in 
L^xyz3s  ^°r  t*ie  P°*nt  groups.)  This  does  not  automatically 


8 


i 


tell  us  which  strains  are  associated  with  which  representations  but 
does  restrict  the  list  of  possible  representations.  For  a given  point 
group  this  list  may  also  be  obtained  indirectly  if  the  list  for  a "parent" 
group  is  known;  for  example,  the  tetragonal  point  group  (D^)  list  may 
be  obtained  from  the  cubic  point  group  (0^)  list  by  using  the 
compatibility  tables.  An  overall  chart  of  the  subgroup  decomposition 
of  the  32  point  groups  is  shown  in  Appendix  D. 

Using  cubic,  hexagonal,  and  tetragonal  crystals  as  examples,  we  now 
show  how  specific  strains  are  associated  with  specific  lower- symmetry 
point  groups. 

B.  Cubic  Crystals  (0^) . 

For  cubic  crystals  we  may  decompose  e as 


+ d 


10  0 
0 10 
0 0 1 

0 0 o' 
0 0 1 
0 1 0 


+ b 


+ f 


1 0 0 
0 I 0 
0 0 2 

0 0 1 
0 0 0 
10  0 


+ c 


+ h 


1 0 0 
0 T 0 
0 0 0 

0 10 
10  0 
0 0 0 


with  a = (en  ♦ e22  ♦ e33)/3 


b = (2e33  - en  - e22^6 
c = (ell  ‘ e22)/2 


(15a) 


(15b) 


(15b) 


d = e 


23 


f = e 


13 


h = e 


12 


Using  the  procedure  of  Appendix  E,  one  finds  that  the  first  matrix 
in  eq  (15a)  transforms  like  x2+y2+z2,  i.e.,  like  P1  (the  correspondence 
between  functions  and  irreducible  representations  can  be  obtained  from 
standard  tables);  the  second  and  third  matrices  transform  like 
2z2-x2-y2  and  x2-y2,  i.e.,  like  the  two  partners  of  P12;  the  last 
three  matrices  transform  like  yz,  xz,  and  xy,  i.e.,  like  the  three 
partners  of  P 2 . This  is  consistent  with  Table  C-l  for  0^. 

For  actual  computation  it  is  probably  more  convenient  to  decompose 


e as 


a 


'l  0 0 
0 10 
0 0 1 


+ b 


T o o' 
o T o 
0 0 2 


♦ c 


"2  0 o' 
o T o 
0 0 I 


"o  i r 

o T I* 

o i T 

+d 

10  1 

+f 

T o i 

+h 

1 0 I 

1 1 0 

T i o 

II  0 

- J 

L 

(16a) 


with  a 
b 
c 


(ell  + e22  + ®33)/3 
(e33  “ e22)/3 
(ell  - e22)/3 


(16b) 


d = (e12  ♦ e23)/2 

f * (*23  ‘ ®13)/2  (16b) 

h = (e12  - e13)/2 

Using  eq  (13)  one  sees  that  only  the  first  matrix  in  (16a)  is  associated 
with  a change  in  volume;  this  matrix  transforms  like  and  represents 
hydrostatic  (cubic)  strain.  The  second  and  third  matrices  each  transfora 
like  P^2  and  are  associated  with  tetragonal  distortions;  since  the 
third  matrix  is  related  to  the  second  by  simply  interchanging  the  x and 
z axes,  no  new  computation  is  required.  Similarly,  the  last  three 
matrices  in  (16a)  each  transform  like  J1  ' and  are  associated  with 
strains  along  the  long  diagonals  ^11^ , (lll^  , and  (fll^  , respectively, 
i.e.,  trigonal  distortions;  again,  only  one  actual  computation  is 
needed,  the  remaining  two  being  determined  by  interchange  of  appropri- 
ate axes.  Thus,  for  cubic  crystals,  an  arbitrary  strain  may  be 
written  as  a linear  combination  of  hudrostatic  (cubic),  tetragonal, 
and  trigonal  strains. 

Examples  of  specific  distortions  of  cubic  crystals: 

1.  Cubic  — *►  tetragonal  (D^)* 

(a)  General  tetragonal. 
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*o  -y 


Figure  1.  Cubic  — ► general  tetragonal. 


e = e 


10  0 


0 10 


0 0 oc 


(17a) 


with  e aQ  = At  and  DC  = Ah/ At.  This  is  most  conveniently  written  as 


=»  »1  »2 

e = Cj  e1  + c,  e 


(17b) 


1 0 0 


e1  = e 010 


0 0 1 


10  0 

e2  = e 0 I 0 


0 0 2 


(17c) 


giving 


c,  = (o(-2)/3  , c,  = (oC+  l)/3. 


(17d) 


Thus,  for  an  arbitrary  Ah/ At  ratio  a component  is  required  and 
the  star-of-lt  theorem  does  not  hold;  there  are  two  "operator  functions". 
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Cj  (x2+y2+z2)  and  c2  (2z2-x2-y2)  transforming  like  and  P12, 
respectively. 

(b)  " P12"  tetragonal. 

When  Ah  is  exactly  2 At  (i.e.,  «C  = 2)  we  see  from 


component  is  required  so  that  volume 


is  preserved  and  the  star-of-k  theorem  holds;  the  operator  function 
is  then  2z2-x2-y2  which  transforms  like  fl-. 


2.  Cubic  ->■  orthorhombic  (D„,  ) 


"Even" 


Figure  2.  Cubic  — *•  orthorhombic  "even1 


with  e aQ  = At.  The  operator  function  is  x2-y2,  second  partner  of  P 
There  is  no  P^  component  (volume  is  preserved);  thus,  the  star-of-ic 
theorem  holds. 

Note:  We  may  represent  this  strain  as 


1 0 0 


2 0 0 


e = c^e  010 
0 0 2 


0 10 

0 0 I 


with  Cj  = 1/3,  = 2/3  so  that,  if  one  already  has  a computer  program 

for  the  distortion  associated  with  the  first  matrix  of  eq  (19),  no 
new  programming  is  required. 

(b)  "Uneven". 


-AS 

-/ 


A-  -y 


-/SAt 


Figure  3.  Cubic  —>■  orthorhombic  "uneven". 


10  0 


e = e O/3O 

0 0 oc 


with  e aQ  = 4*>  °<e  aQ  = As.  We  may  rewrite  eq  (20)  as 


1 0 0 


10  0 


e = c,e  10  10  1+  c e 010  +c  e 010  (211 

2 & 


2 0 0 


with  c 


general  orthorhombic  "uneven"  case  c^O,  i.e.,  a component  is 
contained  in  the  strain  and  the  star-of-k  theorem  does  not  hold. 


Note:  If  As  * (1-/3)  At,  i.e.,  cC=ly3  so  that  <^=0,  no  Pj 
component  is  contained  (volume  is  preserved)  and  the  star-of-k  theorem 


holds.  For  this  case  c?»(2/3-l)/3,  c*=(l+/5)/3 


3.  Cubic  — trigonal  (D_ ,) 


Figure  4.  Cubic  —►trigonal 


with  2 eV?  aQ  = A,t.  The  operator  function  is  yz+xz+xy,  the  sum  of 
the  three  partners  of  T ' . There  is  no  P component  (volume  is 


preserved);  thus,  the  star-of-k  theorem  holds 


I 

I 


I 

I { 


I 


C.  Hexagonal  Crystals  (D,_  ) 

oh 

Axes  are  defined  in  Fig.  5.  This  is  consistent  with  Koster  et  al.2 


A ^ 

Figure  5.  Coordinates  for  hexagonal  crystals. 

For  hexagonal  crystals  it  it  probably  most  convenient  to  write  an 
arbitrary  strain  as 


e”  = a 

1 0 0 
0 10 

+ b 

0 0 0 
0 0 0 

+ c 

10  0 
0 I 0 

+ d 

0 0 0 

0 10 
10  0 

+ f 

0 0 1 

0 0 1 
0 0 0 

+ h 

0 0 0 

0 0 0 
0 0 1 

0 0 0 

1 0 0 

0 1 0 

^"Properties  of  the  Thirty-Two  Point  Groups",  G.F.  Koster,  J.O.  Dimmock, 
R.G.  Wheeler,  and  H.  Statz  (M. I.T.  Press,  1963). 
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with  a = (e^  + ®22)/2 


b = e 


33 


C = Ceu  - e22)/2 


C23b) 


d = e 


12 


f = e 


13 


h = e 


23 


The  matrices  in  eq  (23a)  transform  like  x2+y2(  Pj),  z2(  P^, 
x2-y2(P6),  xy(P6),  xz(r5),  and  yz(P5);  this  is  consistent  with 
Table  C-l  for  D^.  The  first  two  strains  in  (23a)  preserve  the 

hexagonal  symmetry;  the  third  strain  is  associated  with  an  orthorhombic 

distortion;  the  last  three  strains  are  associated  with  monoclinic 
distortions  (the  xy  distortion  forms  a monoclinic  P lattice,  while  the 
xz  and  yz  distortions  form  monoclinic  C lattices). 

Note:  A strain  for  which  b = -2a  in  eq  (23a)  will  be  volume- 
preserving (see  eq  (13))  but  will  still  contain  so  that  the  star- 

of-k  theorem  will  not  hold. 

Examples  of  specific  distortions  of  hexagonal  crystals: 

1 . Hexagonal  — ► orthorhombic  (D2h) • 


17 


Figure  6.  Hexagonal orthorhombic. 


10  0 


e = e 0 1 0 


0 0 0 


with  e aQ  = At.  The  operator  function  is  x^-y^,  a partner  of  f^. 
Since  there  is  no  component  the  star-of-k  theorem  holds.  In 
Fig.  6 the  points  1,2, 3, 4 form  the  base  of  the  orthorhombic  C cell. 
Before  distortion  the  sides  of  this  cell  have  the  ratio  'P  required 
for  a hexagonal  crystal;  after  distortion  the  ratio  becomes 
VP  ao(l-e)/aQ(l+e)  = V?(l-2e). 


a ' = a , a ' = a (1+  fPe/2)  (25b) 

1 o 2 o 

The  operator  function  is  xy,  a partner  of  P . Since  there  is  no  P 

6 1 

component  the  star-of-k  theorem  holds.  In  Fig.  7 the  dotted  lines 
connecting  points  1,2, 3, 4 form  the  base  of  a monoclinic  P crystal  with 
angle  1,4,3  = /3  • The  vertical  sides  of  the  monoclinic  P are  rectangles 
with  height  cq  (unperturbed  hexagonal  height) . 


3.  Hexagonal  -*  monoclinic  C (^2^)  - xz. 

In  Fig.  6 the  points  1,2, 3, 4 form  the  base  of  an  orthorhombic  C 
cell.  Under  xz  distortion  the  vertical  (z  direction)  face  with  bottom 
side  1,2  is  distorted  from  a rectangle  to  a parallelogram  and  the  solid 
with  base  1,2, 5, 6 and  vertical  sides  forms  a primitive  cell  for  the 
monoclinic  C cell  (with  base  1,2, 3, 4).  Details  of  the  geometry  are 
given  in  Appendix  F,  paragraph  3. 


e = e 


0 0 1 
0 0 0 
1 0 0 


(26) 


The  operator  function  is  xz,  a partner  of  Pg.  Since  no  P^  component 

— fc 

is  contained  the  star-of-k  theorem  holds. 


4.  Hexagonal  -*•  monoclinic  C (C^)  - yz. 

Again  referring  to  Fig. 6 the  points  1,2, 3, 4 form  the  base  of  an 
orthorhombic  C cell.  Under  yz  distortion  the  vertical  (z  direction) 
face  with  bottom  side  2,3  is  distorted  from  a rectangle  to  a parallelogram. 
As  in  the  previous  case,  the  solid  with  base  1,2, 5, 6 and  vertical  sides 


20 


forms  a primitive  cell  for  the  monoclinic  C cell  (with  base  1,2, 3, 4). 
Details  of  the  geometry  are  given  in  Appendix  F,  paragraph  4. 


0 0 0 
0 0 1 
0 1 0 


(27) 


The  operator  function  is  yz,  a partner  of  Since  no  I"1  component 

is  contained  the  star-of-k  theorem  holds. 

Note:  the  primitive  cells  for  the  xz  and  yz  distortions  do  not 

have  the  same  shape. 


D.  Tetragonal  Crystals  (D^). 

| 

For  tetragonal  crystals  (with  z-axis  unique)  an  arbitrary  strain  can 
be  written  as 


"e  = a 

"l  o o' 
0 10 

♦ b 

'o  0 o' 
0 0 0 

+ c 

1 0 o' 
o T o 

0 0 0 

_0  0 1^ 

0 0 0^ 

(28a) 


0 1 0 

0 0 1 

0 0 0 

1 0 0 

+ f 

0 0 0 

+ h 

0 0 1 

0 0 0, 

1 0 0 

0 10 

with  a = (e..  + el/2 


(28b) 


d = e. 


f = e. 


h = e 

23 

The  first  matrix  in  eq  (28a)  transforms  like  x2  + y2  ( pp,  the  second 

as  z2  ( pj)»  the  third  as  x2  - y2  ( J^),  the  fourth  as  xy  ( T^),  and 

the  last  two  as  xz  and  yz,  respectively  (the  two  partners  of  P5) ; this 

is  consistent  with  Table  C-l  for  D,,  . 

4h 

Note:  A strain  for  which  b * -la  in  eq  (28a)  will  be  volume- 

preserving (see  eq  (13))  but  will  stilt  contain  so  that  the  sts. 
of-k  theorem  will  not  hol.l. 

Examples  of  specific  distortions  of  tetragonal  crystals: 

1.  Tetragonal — * orthorhombic  u>_.). 

2h 


(a)  "Even". 


Te »| 


-At 


Figure  8.  Tetragonal —►orthorhombic  "even". 


10  0 


0 0 0 


with  et0  = ^ ' 1116  °Perator  function  is  x2  - y2  ( P2).  Since  there 

is  no  Pj  component  the  star-of-k  theorem  holds.  Volume  is  preserved 
(b)  "Uneven". 


Figure  9.  Tetragonal— ► orthorhombic  "uneven" 


(1  VSJ/2 


If  As  = (1-/3)  (s0/tQ)At;  i.e.,  oC=  (1-/3),  we  see  from  either  (30) 
or  (31)  that  volume  is  preserved.  However,  there  are  still  P^ 
components  (unless  oC  = 0 and  1;  i.e.,  the  orthorhombic  "even"  case) 


so  that  the  star-of-k  theorem  does  not  hold.  The  three  operator 


functions  associated  with  eq  (31)  are  x2+y2,  z2,  and  x2-y  , respectively 
or  ct  P^l)  + c2  Pj(2)  ♦ C3  p2. 

2.  Tetragonal  — * orthorhombic  C (^2^)- 


Figure  10.  Tetragonal  — y orthorhombic  C. 


e = e 


0 10 
10  0 
0 0 0 


(32) 


The  operator  function  is  xy  ( P^).  Since  no  P^  component  is  contained 

-X 

the  star-of-k  theorem  holds.  The  solid  with  base  1,2, 3, 4 and  vertical 
(z-direction)  sides  is  a primitive  cell  of  the  orthorhombic  C cell  with 
base  1,5, 7, 3 and  vertical  sides.  All  vertical  sides  are  s0  (undeformed 


tetragonal  height).  The  new  primitive  cell  has  sides  t_,  t , s (to 


first  order);  angle  1,2,3  is  90°-2e.  The  orthorhombic  C cell  has  base 
sides  1,3  (V?t0(l-e))  and  3,7  flfTt  (1+e)) . 


3.  Tetragonal  monoclinic  P (C_.  ) - xz. 

2h 


/ 

/ 


Figure  11.  Tetragonal  monoclinic  P. 

I 0 0 l| 


e = e 0 0 0 


1 1 0 0| 

The  operator  function  is  xz,  a partner  of  I1  . Since  no  component 
is  contained  the  star-of-k  theorem  holds.  After  distortion  the  solid 
with  base  1,2, 3, 4 and  vertical  sides  (y  direction)  is  a monoclinic  P 
cell  with  side  1,4  * sQ,  side  1,2  ■ tQ,  and  vertical  side  tQ;  angle 
2,3,4  is  90°-2e. 

Note  1:  The  mcnoclinic  P cell  used  here  is  "special"  in  that  two 

sides  have  the  same  length  tQ.  No  new  Bravais  lattice  is  thereby 
created,  however,  since  this  cell  has  no  additional  symmetry  elements 


i 


over  and  above  those  associated  with  a cell  of  identical  base  but 


vertical  height  not  equal  to  t0. 

Note  2:  From  Fig.  11  one  sees  why  there  is  no  separately  named 

monoclinic  C Bravais  lattice  with  face-center  atoms  in  the  non-90°- 
angle  face  (base  2, 4, 6, 5);  such  a lattice  can  always  be  considered  as 
monoclinic  P with  base  1,2, 3, 4. 

Note  3:  The  tetragonal  —*■  monoclinic  P ^^3  (yz)  case  is  identical 

to  the  xz  case  with  the  roles  of  the  x and  y axes  interchanged. 


I 

j 

I 


t 
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APPENDIX  A 


EQUIVALENT  LABELS  FOR  IRREDUCIBLE  REPRESENTATIONS 


TABLE  A-l. 

EQUIVALENT  LABELS  FOR  THE 

CUBIC  GROUP 

BSWa 

Koster  et  al.b 

Tinkham0 

r, 

Ti* 

A1 

r2 

P2+ 

A2 

^12 

E 

r!5' 

V 

T1 

r2S' 

Ts* 

T2 

<v 

<V 

r2- 

^12/ 

V 

ris 

r4- 

^25 

r5- 

a Ref  1,  Table  I.  Koster  (ref  13,  Table  XXXI)  and  Slater  (ref  14, 
Table  A3-20)  also  use  the  BSW  notation. 

b Ref  2,  Table  87.  This  is  also  known  as  the  Bethe  notation. 

c Ref  12,  Table  for  0,  App  B,  p.  329.  Tinkham  gives  labels  for  0 
(proper  rotations  only).  This  is  also  known  as  the  Mulliken  notation. 

1.  L.P.  Bouckaert,  R.  Smoluchowski,  and  E.  Wigner,  Phys.  Rev.  50^ 

58  (1936). 

2.  "Properties  of  the  Thirty-Two  Point  Groups,"  G.F.  Koster, 

J.O.  Dimmock,  R.G.  Wheeler,  and  H.  Statz  (M.I.T.  Press,  1963). 

12.  "Group  Theory  and  Quantum  Mechanics,"  M.  Tinkham  (McGraw-Hill, 

1964) . 

13.  G.F.  Koster,  "Space  Groups  and  Their  Representations"  in  "Solid 
State  Physics,  Vol.  5,"  ed.  F.  Seitz  and  D.  Turnbull  (Academic 
Press,  1957) . 

14.  "Quantum  Theory  of  Molecules  and  Solids,  Vol.  2,"  J.C.  Slater 
(McGraw-Hill,  1965). 
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TABLE  A-2 . EQUIVALENT  LABELS  FOR  THE  HEXAGONAL  GROUP  (D6h) 

Koster  et  al.a 

Slaterb 

Tinkham0 

Ti<- 

A1 

iv 

r2  v 

A2 

h f 

r4v 

B1 

r4<- 

B2 

r5<- 

E1 

r6+ 

r5* 

E2 

Ti- 

V 

r2- 

r2- 

r3- 

r4- 

Pr 

r3- 

r5- 

r6- 

r6- 

r5- 

aRef  2,  Table  12.  Koster  (ref 
(unlabeled)  in  the  same  order. 

bRef  14,  Table  A3-2. 

13,  Table  XXVI) 

lists  the  representations 

cRef  12,  Table  for  Dg,  App  B, 

p.  327.  Tinkham 

gives  labels  for 

(proper  rotations  of  D^) . 

2.  "Properties  of  the  Thirty-Two  Point  Groups,"  G.F.  Koster, 

J.O.  Dinunock,  R.G.  Wheeler,  and  H.  Statz  (M.I.T.  Press,  1963). 


12.  "Group  Theory  and  Quantum  Mechanics,"  M.  Tinkham  (McGraw-Hill, 
1964) . 

13.  G.F.  Koster,  "Space  Groups  and  Their  Representations"  in  "Solid 
State  Physics,  Vol.  5,"  ed.  F.  Seitz  and  D.  Turnbull  (Academic 
Press,  1957). 

14.  "Quantum  Theory  of  Molecules  and  Solids,  Vol.  2,"  J.C.  Slater 
(McGraw-Hill,  1965). 
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TABLE  A-3.  EQUIVALENT  LABELS  FOR  THE  TETRAGONAL  GROUP  (D4h) 
This  worka  Koster  et  al.b  Kosterc  Slaterd  Tinkhame 

P1  ?1*  M1  A. 


a rt  here  is  X^(M^)  of  BSW,  ref  1,  Table  V.  This  notation  is  also 
used  in  ref  8 and  9. 


D Ref  2,  Table  40. 

c Ref  13,  Table  XIV. 

d Ref  14,  Table  A3-74. 

e Ref  12,  Table  for  D^,  App  B,  p.  327.  Tinkham  gives  labels  for  D 
(proper  rotations  or  D^) . 4 

1.  L.P.  Bouckaert,  R.  Smoluchowski , and  E.  Wigner,  Phys.  Rev.  50, 

58  (1936). 

2.  "Properties  of  the  Thirty-Two  Point  Groups,"  G.F.  Koster, 
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9.  D.M.  Gray  and  A.M.  Gray,  Phys.  Rev.  B,  J_4,  669  (1976). 

12.  "Group  Theory  and  Quantum  Mechanics,"  M.  Tinkham  (McGraw-Hill, 
1964). 

13.  G.F.  Koster,  "Space  Groups  and  Their  Representations"  in  "Solid 
State  Physics,  Vol.  5,"  ed.  F.  Seitz  and  D.  Turnbull  (Academic 
Press,  1957). 

14.  "Quantum  Theory  of  Molecules  and  Solids,  Vol.  2,"  J.C.  Slater 
(McGraw-Hill,  1965). 
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APPENDIX  B 

PROOF  OF  COROLLARIES  TO  THE  "STAR-OF-k"  THEOREM 

Corollary  1:  For  cubic  crystals  AVolume=0  is  the  necessary  (N) 

and  sufficient  (S)  condition  for  AE_  (etc.)=0. 

r 

N:  From  the  proof9  of  the  main  theorem,  the  absence  of  a 

component  in  e”  is  the  N and  S condition  that  AEp=0.  The  absence  of 
rx  implies  that  the  sum  of  all  the  operations  of  the  cubic  group  applied 
to  * equals  zero9.  Since  these  different  rotations  must  all  produce  the 
same  volume  change,  AVolume=0.  Thus  AE^=0  =$  no  P,  Avolume=0; 

r 1 

then  AVolume=0  is  necessary  for  /±E  =0.  Q.E.D. 

F * 

S:  For  a cubic  crystal  a P^  strain  (a  strain  which  preserves  the 

cubic  symmetry)  must  always  have  a volume  change  associated  with  it;  thus 
AVolume  ^fcO.  Therefore,  AVolume=0  =^no  P^,  and,  as  no 
pj  — v-  AEp=0,  AVolume=0  is  sufficient  for  AEp=0.  Q.E.D. 

Corollary  2:  For  non- cubic  crystals  AVolume=0  is  N but  not  S for 

AEp(etc.)=0. 

N:  Proof  is  exactly  as  in  corollary  1. 

S:  For  non-cubic  crystals  one  can  preserve  the  symmetry  without 

changing  volume.  Consider  a tetragonal  crystal  txtxs,  for  example. 

The  distortion  t-*t(l+a),  s -*s(l-2a)  preserves  volume  to  first  order. 

The  tetragonal  symmetry  is  preserved;  one  has  merely  changed  the  s/t  ratio. 
Thus  Ayolume=0  is  not  sufficient  for  AEp=0.  Q.E.D. 

^D.M.  Gray  and  A.M.  Gray,  Phys.  Rev.  B,  14_,  669  (1976). 
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APPENDIX  C 


USE  OF  [l^  ls  TO  GENERATE  THE  REPRESENTATIONS  CONTAINED  IN  e 

Definition:  Given  two  representations  P and  p of  the  same  group 

m n 

.H  with  basis  functions  and  respectively,  the 

direct-product  matrices  D(R)  = Dm(R)xDn(R),  with  R any  operator  in 
% , form  a representation  of M with  basis  functions  . . If  the 

two  sets  of  functions  0,  J*  belong  to  the  same  representation,  Pn  of 
Jfcf.  then  the  "symmetrized- square"  representation  frjs  will  be 
associated  with  functions  of  types  ^ and  h * A r the 

"antisymmetrized- square"  representation  cr„]  _ will  be  associated  with 

*•  a 

functions  of  type  01  ?!  - tj  K 
Traces: 

Direct-product  XmXn(R)  = Xm(P)  xX"(R) 

Symmetrized-square  2 X^n*n^s(R)  = CXn(R)]2  ♦ 

Antisymmetrized- square  2 X^nXn^a(R)  = fXn(R)l2  - ^Cn(R2) 

An  arbitrary  strain  can  be  represented  by  six  independent  e00* 
as  indicated  in  eq  (14)  of  the  main  text,  thus  six  linearly  independent 
functions  will  be  needed;  these  may  be  taken  as  x^,  , xy,  xz,  and 

yz.  Since  basis  functions  for  a direct-product  representation  are  the 
various  products  of  the  basis  functions  for  the  original  representations, 
the  original  P ^ with  basis  functions  x,y,  and  z will  generate 
xX,xY,xZ,yX,yY,yZ,zX, zY,  and  zZ.  We  desire  the  symmetric  functions 
xX  (x^) , xY+yX  (xy),  etc.  These  six  functions  will  be  basis  functions 
for  the  representation  comprising  the  symmetrized-square  of  P 
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(xY-yX,  etc.  will  be  associated  with  the  antisymmetrized  square.)  The 
6-dimensional  Q Txyz  1 s will  be  reducible;  the  contained  irreducible 
representations  will  be  just  the  irreducible  representations  required  to 
express  the  six  strain  functions  and  thus  just  the  required  set  to 
express  an  arbitrary  strain.  (See  also  ref  12,  pp  43-47  and  ref  15, 
pp  11-44  to  11-48.) 


TABLE  C-l.  IRREDUCIBLE  REPRESENTATIONS3  CONTAINED  IN  fp  1 

L1  xyz  4 s 


Point  group3 

rxyz 

[ r 1 

Cubic  (0h) 

r!5 

P1  * ri2  * 

^25/ 

Hexagonal  (D^) 

r2-  * r5- 

2 r i * p5  * 

P 

6 

Tetragonal  (D.,  ) 
4h 

r4'  * rs/ 

2 Pj  * r2  * 

r3*  p5 

0 

r4 

r + r + 

1 3 

1 

Td 

r5 

*1  * P3  * 

rs 

Th 

r4- 

Pi  * r2  * 

r3  * r„ 

T 

r4 

Pi  * p2* 

r3*  P. 

*2"Group  Theory  and  Quantum  Mechanics,"  M.  Tinkham  (McGraw-Hill,  1964). 
^Lecture  notes  on  "Theory  of  the  Properties  of  Ionic  Complexes  in 
Crystals,"  Frank  Ham,  RPI,  1969  (unpublished). 


33 


v i 

l 


APPENDIX  D 


SUBGROUP  DECOMPOSITION  OF  THE  32  POINT  GROUPS 


Figure  12 


Subgroup  decomposition  of  the  32  point  groups.  A heavy- 
line  indicates  that  the  subgroup  is  not  invariant. 


In  Fig.  12,  all  of  the  32  point  groups  are  illustrated  along  with 


the  possible  subgroups  of  each  of  these  groups.  Light  lines  indicate 


invariant  subgroups,  heavy  lines  indicate  non-invariant  subgroups. 

(Fig.  12  is  taken  directly  from  ref  2;  compatibility  tables  connecting 
these  groups  are  given  therein.) 


"Properties  of  the  Thirty-Two  Point  Groups",  G.F.  Koster,  J.O.  Dimmock 
R.G.  Wheeler,  and  H.  Statz  (M.I.T.  Press,  1963). 


APPENDIX  E 

TRANSFORMATION  PROPERTIES  OF  OPERATORS  DETERMINED  BY  TRANSFORMATION 
PROPERTIES  OF  BASIS  FUNCTIONS 

When  the  desired  strain  e decomposes  into  ^ c^  ^ e1*  it  is 
not  always  apparent  how  the  various  e^4*  transform.  If  one  knows  how 
basis  functions  for  various  representations  transform  it  may  be  easier 
to  first  associate  functions  f*  (r)  with  the  various  el0t  . We  now  show 


that  for 


fj^  (r)  3 rT  eixr. 


(E-l) 


then  the  eio<  transform  like  the  f*  (r) . Here  "r  is  the  column  vector 


fi(i 


and  r^  is  the  row  vector 

rT  = (x,y, z) . 

Proof:  The  necessary  and  sufficient  condition  that  a matrix  e 

transform  like  representation  i is  that 

PR  (eioc)  PR-|  =JADXn^®iX  (E'2) 
for  all  operators  R of  the  pertinent  group.  Multiply  eq  (E-2)  by  rT 


from  the  left  and  by  r from  the  right: 


r ' e1**'  - IAD^  (R)  f^  (r),  (E-3) 

using  (E-l)  and  PR-l  r = r'.  From  the  definition  of  f ^ (r)  in  (E-l), 
the  left  side  of  (E-3)  is  just  f^(r');  thus  (E-3)  becomes 
fV(r')  = f£  (R-1r)  -^D^R)  f^(r). 

If  this  is  repeated  for  of,  = 1,2. . . ,d  where  d is  the  dimensionality  of 
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representation  i,  then 


fj  (R-l?) 


f|  (R“*r) 


D1(R) 


ffJ(r) 


(E-4) 


Since  eq  (E-4)  is  true  for  all  R,  this  is  just  the  condition  that  the 
f*  (r)  transform  like  representation  i,  i.e.,  that  the  f*  (r)  defined 
by  eq  (E-l)  are  basis  functions  for  representation  i.  Thus,  the 
transformation  properties  of  the  e*-*  are  the  same  as  those  of  the 
fi  defined  by  eq  (E-l),  Q.E.D. 


< 


APPENDIX  F 

HEXAGONAL  CASE:  GEOMETRIC  DETAILS 

1.  ("^  strains. 

(a)  For  a strain  associated  with  x2+y2  (en  = e22  = e,  all  other 
eij  = every  vector  in  the  xy  plane  will  simply  be  extended  by  a 
factor  (1+e).  Thus,  the  hexagonal  symmetry  is  maintained;  one  merely 
decreases  the  c/a  ratio  (for  positive  e). 

(b)  For  a strain  associated  with  z2  (e33  = e , all  other  e^  = 0) 
the  hexagonal  symmetry  is  clearly  maintained;  one  merely  increases  the 
c/a  ratio  (for  positive  e) . 

2.  The  strain  xy. 

For  a strain  associated  with  xy,  e12  = e21  = e,  all  other  e^  = 0. 
Referring  to  Fig.  7 in  the  main  text 

a^Q  = ai,  a*20  = (a/2)i  + (/?a/2)*. 

Using  eq  (12)  of  the  main  text  these  become 
a^  = ai  + eaj 

t2  = a (1/2  + V?e/  2)i  + a(f?/2  + e/2)j\ 

To  first  order  a^=a,  a2=a  (l+i3e/2),  and  the  angle  between  a^  and  i*2 
is  60°-3e/2.  The  detailed  geometry  of  this  distortion  is  indicated 
in  Fig.  13. 
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3.  The  r*s  strain  xz. 

In  Fig.  6 of  the  main  text  take  the  origin  at  point  2,  x-axis 
toward  point  1,  y-axis  toward  point  3,  and  z-axis  vertically  up  from 
the  plane  of  the  figure.  Before  distortion  the  points  1,2, 3, 4 form  the 
base  of  an  orthorhombic  C solid.  Under  xz  distortion  a non-90°  angle 
is  produced  in  the  vertical  side  with  bottom  edge  2,1.  Define  as 
the  undistorted  vector  from  the  origin  to  point  1,  a^^  as  the  vector 
from  the  origin  to  point  3,  and  c^  as  the  pertinent  vertical  vector. 
Under  xz  distortion 

•4  A A 

al  = V + eaok 
*3  = *30  = 

U A A 

c = ecQi  ♦ cQk 

so  that,  to  first  order,  a2=a0,  a^V^kg,  c=cQ,  and  the  angle /3  between 
a’j  and  c is  90°-2e.  A perspective  drawing  is  shown  in  Fig.  14.  In 


I 


■ 
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monoclinic  C cell  is  outlined  by  points  1-4,  7-10;  ot*)f=900t  /3  =90°-2e. 
The  primitive  cell  is  outlined  by  points  4,6,3,11,  etc,  side  4,6  = aQ, 
side  6,3  = aQ,  side  3,9  = c0,  angle  6,3,11  = 120°,  angle  6,3,9  = 90°-e. 

4.  The  P5  strain  yz. 

In  Figure  6 of  the  main  text  take  the  origin  at  point  2,  x-axis 
toward  point  1,  y-axis  toward  point  3,  and  z-axis  vertically  up  from  the 
plane  of  the  figure  (exactly  as  in  the  xz  case) . Before  distortion 
the  points  1,2, 3, 4 form  the  base  of  an  orthorhombic  C solid.  Under  yz 
distortion  a non-90°  angle  is  produced  in  the  vertical  side  with  bottom 
edge  2,3.  Define  a1Q,  a^,  and  cQ  as  in  the  previous  section.  Under 
yz  distortion 

-k  A A 

ai  = aio  = V 

^3  = + 

A A 

c = ec  i + c k 

o o 

so  that,  to  first  order,  a^=aQ,  a.^^a.Q,  c=cQ,  and  the  angle /3  between 
ifj  and  "c  is  90° -2e.  A perspective  drawing  is  shown  in  Fig.  15. 
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monoclinic  C cell  is  outlined  by  points  1-4,  7-10;  Ct  = y * 90°, 
yj  =90°-2e.  The  primitive  cell  is  outlined  by  points  4,6,1,11  etc 
side  1,6  = aQ,  side  1,11  = aQ,  side  1,7  = cQ,  angle  4,6,1  = 120°. 

angle  6,1,7  = 90°-\l?e. 
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